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Graphene has raised high expectations as a low-loss plasmonic material in which the plasmon 
properties can be controlled via electrostatic doping. Here, we analyze realistic configurations, 
which produce inhomogeneous doping, in contrast to what has been so far assumed in the study 
of plasmons in nanostructured graphene. Specifically, we investigate backgated ribbons, co-planar 
ribbon pairs placed at opposite potentials, and individual ribbons subject to a uniform electric field. 
Plasmons in backgated ribbons and ribbon pairs are similar to those of uniformly doped ribbons, 
provided the Fermi energy is appropriately scaled to compensate for finite-size effects such as the 
divergence of the carrier density at the edges. In contrast, the plasmons of a ribbon exposed to a 
uniform field exhibit distinct dispersion and spatial profiles that considerably differ from uniformly 
doped ribbons. Our results provide a road map to understand graphene plasmons under realistic 
electrostatic doping conditions. 



Plasmons^ -the collective oscillations of conduction 
electrons in metals- are capable of confining electro- 
magnetic energy down to deep sub- wavelength regions. 
They can also enhance the intensity of an incident light 
wave by several orders of magnitude. These phenom- 
ena are the main reason why the field of plasmonics is 
finding a wide range of applications that include single- 
molecule sensing, 2 nonlinear optics,^ and optical trapping 
of nanometer-sized objects P 

Recently, confined plasmons have been observed and 
spatially mapped in doped grapheneP The level of dop- 
ing in this material can be adjusted by exposing it to 
the electric fields produced by neighboring gates. Elec- 
trostatic doping has actually been used to demonstrate 
plasmon-frequency tunability^ and induced optical mod- 
ulations in the THz^and infrarecP response of graphene. 

The two-dimensional (2D) band structure of pristine 
graphene consists of two cones filled with valence elec- 
trons and two empty inverted cones joining the former at 
the so-called Dirac points, which mark the Fermi level. 
Extra electrons or holes added to this structure form a 2D 
electron or hole gas that can sustain surface plasmons 
Compared to noble-metal plasmons, graphene modes are 
believed to be long-lived excitations. 10 But most impor- 
tantly, their frequency can be controlled via the above- 
mentioned electrostatic doping. 5 7 For example, in homo- 
geneous suspended graphene, a perpendicular DC elec- 
tric field £ applied to one side of the carbon sheet is 
completely screened by an induced surface charge den- 
sity —en = £/47r, and this in turn situates the Fermi 
level at an energy Ep = hvp y^^\ n \ relative to the Dirac 
points. 11 Here, vp ~ 10 6 m/s is the Fermi velocity of 
graphene. 

Plasmons in doped graphene nanostructures have been 
generally studi ed by assuming a uniform doping electron 
density n ! 12 | 1S | But in practice n is inhomogeneous and 
depends on the actual geometrical configuration. For 
example, in a ribbon of width W placed at a distance 
d from a planar biased backgate, n shows a dramatic 
pileup near the edges, as it is well known in microstrip 
technologypH^ leading to divergent local Ep levels, as 



shown in the inset of Fig.[TJa) for various values of W/d. 
More precisely, n diverges as oc 1/y/x with the distance 
to the edge x (and hence, Ep oc l/x 1//4 ). The average of 
the Fermi energy over the ribbon area [(|i^p|), see Fig. 
[TJa)] is very different from the W/d > 1 limit (Ef = 
hv F yJ\V\/4ed) and diverges as 18 y/dJW / \og{2d/W) in 
the narrow ribbon limit (W <C d) for constant bias po- 
tential V . Therefore, the question arises, how different 
are the plasmon energies and field distributions in ac- 
tual doped graphene nanostructures compared to those 
obtained for uniformly doped graphene? 

Here, we analyze plasmons in doped graphene ribbons 
under different geometrical configurations. Specifically, 
we study backgated single ribbons, co-planar parallel rib- 
bon pairs of opposite polarity, and single ribbons im- 
mersed in a uniform external electric field. For sim- 
plicity, we describe the frequency-dependent conductiv- 
ity of doped graphene in the Drude model as a(uj) = 
(ie 2 \Ep\/7rh 2 )/(uj + ry), where 7 <C u is a small relax- 
ation rate. The doping electron density n is obtained 
from electrostatic boundary-element calculations, while 
the plasmon frequencies are computed using a discrete 
surface-dipole approximation (DSDA), as explained in 
Appendix [Cj 

Backgated ribbons.- A first conclusion extracted from 
Fig. [TJa) is that the level of doping, quantified in the 
Fermi energy, is not well described by the simple capaci- 
tor analysis of the W ^> d geometry. Normalizing to the 
average value (\Ep\), we find Fermi-energy profiles that 
vary between a well of sharp corners for large W/d and 
a smoother, converged shape for small W/d [upper inset 
of Fig. 1(a)]. The former limit corresponds to the ribbon 
in close proximity to the backgate, in which Ep is nearly 
uniform. In contrast, at large separations (d ^> W) we 
find a profile determined by the interaction with a dis- 
tant image, which converges to a well-defined shape up 
to an overall factor (\Ep\) evolving as shown in the main 
plot of Fig. 1(a). 

It is convenient to normalize the ribbon plasmon fre- 
quencies to ujo = (e/h)^/(\Ep\)/W^ so that uj/ujq is a 
dimensionless number, independent of the specific rib- 
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FIG. 1: Electrostatic doping and plasmon modes in backgated 
graphene ribbons, (a) Average Fermi energy as a func- 

tion of width-to-distance ratio W/d, normalized to the value 
obtained in the W ^> d limit. The upper inset shows the 
Ep distribution across the ribbon, normalized to (\Ep\). The 
lower inset shows a sketch of the geometry, (b) Frequency 
uj of the dip olar and q uadrupolar plasmons, normalized to 
cjo — ( e /fr)V (\Ef\) / W , as obtained from the Drude model. 
The insets show the surface charge- density oscillating at fre- 
quency uj and corresponding to these plasmons (vertical axis) 
as a function of position across the ribbon (horizontal axis). 
The dashed curves indicate the W ^> d limit. 



bon width W and gate voltage (i.e., {\Ep\)), as proved 
in Appendix [B] For example, with W = 100 nm and 
(\Ef\) = 0.5 eV, we find huoo = 0.085 eV and a dipole 
plasmon energy huj ~ 0.17eV (wavelength ~ 7.3 jam). 
With this normalization, uj/ujq shows just a mild depen- 
dence on W/d for the dipolar and quadrupolar modes 
[Fig. [IJb)]. The corresponding induced densities (insets) 
are only slightly affected by the change in doping pro- 
file relative to uniform doping (i.e., the average level of 
doping is a dominant parameter, and the effect of edge 
divergences is only marginal). In conclusion, the plasmon 
frequencies and induced densities can be approximately 
described by assuming a uniform Fermi energy in back- 
gated ribbons, thus supporti ng the v alidity of previous 
analyses for this configurationp^^^ although the Fermi 
energy has to be appropriately scaled as shown in Fig. 
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FIG. 2: Plasmons in pairs of co-planar parallel graphene 
ribbons of opposite polarity, (a) Fermi energy distribution 
across pairs of ribbons for different ratios of the ribbon width 
W to the gap distance d. The ribbons are placed at poten- 
tials — V and V, respectively (see inset), and the Fermi energy 
Ef is normalized to the value Ep° — hvF\/\V\/4ed. (b) Fre- 
quency uj of the dipolar and quadrupolar modes as a function 
of W/d, normalized to uj = (e/h)y/(\E F \)/W. Solid (dashed) 
curves correspond to inhomogeneous (uniform) doping. The 
insets show the plasmon charge-density associated with both 
modes (vertical axis) as a function of position across the rib- 
bon on the right (horizontal axis, with the position of the gap 
indicated by an arrow) for W/d = 0.2,1,3, and 10 (curves 
evolving in the direction of the arrows). 



[TJa) to compensate for the effect of finite W/d ratios. 

Two co -planar parallel graphene ribbons.- Two neigh- 
boring ribbons can act both as plasmonic structures and 
as gates. We explore this possibility in Fig. [2j where 
the ribbons are taken to be oppositely polarized. This 
produces doping profiles as shown in Fig. |2|a), which 
evolve from a shape similar to the one obtained for the 
single ribbon of Fig. [l] in the small-ribbon limit at large 
ribbon-pair separations, towards a converged profile near 
the gap in the W/d — >• oo limit. Again, plasmons in this 
structure are very similar to those of neighboring uni- 
formly doped ribbons [see Fig. ^[b)], provided one com- 
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FIG. 3: Plasmons in individual ribbons subject to a uniform 
external electric field £. (a) Fermi energy distribution nor- 
malized to (\Ep\) — 0.6 Hvf y/\S\Je. The inset shows a sketch 
of the geometry, (b) Surface charge associated with the dipo- 
lar plasmon mode (solid curve). The plasmon frequency is 
uj ~ 1.45 cjo, where ujq — (e/K)yJ {\Ef\)/W. The dashed curve 
shows the charge density profile for a uniform doping density 
(i.e., Ep — (\Ep\))- (c) Plasmon dispersion diagram rep- 
resenting the dependence of the density of optical states on 
frequency uj and wave vector parallel to the ribbon k\\. (d) 
Same as (c) for uniform doping. 



pares for the same value of the average Fermi energy 
(\Ef\). The separation dependence of (\Ef\) is shown in 
Appendix [A] Incidentally, plasmons in pairs of uniform 
ribbons have been thoroughly described and the evolu- 
tion of the plasmon frequency with distance explained 
in a recent publication, 20 including the redshift with de- 
creasing d. 



Doping through a uniform electric field.- A disadvan- 
tage of the above doping schemes is the fabrication pro- 
cess involved in adding contacts that allow electrically 
charging the graphene, which can be a source of defects 
in the carbon layer. This could be avoided by doping 
through an external electric field £ produced by either 
distant gates or low-frequency radiation. The graphene 
would then remain globally neutral. This possibility is 
analyzed in Fig. |3j where we consider a ribbon subject 
to a uniform electric field directed across its width. The 
doping profile [Fig. [3^a)] is again exhibiting a divergence 
of Ep at the edges, and it vanishes at the center of the 
ribbon, whe re the doping density changes sign (see Ap- 
pendix _ A3). Following the methods described in Ap- 
we find the relation {\Ep\) = 0.6 hvp^\£\/e. 



A3 



ing dipolar plasmon [Fig. [3^b), solid curve] 



pendix 
The result!: 

displays a large concentration of induced charges near 
the center of the ribbon, in contrast to the dipole plas- 
mon obtained for a homogeneous doping density (dashed 
curve). This inhomogeneous dipole-charge concentration 
is induced by the vanishing of the doping charge density, 
which can be understood as a thinning of the effective 
graphene- layer thickness, similar to what happens near 
the junction of two barely touching metallic structures 
(e.g., two spheres 2 ^). 

We have so far discussed plasmons that are invariant 
along the length of the ribbon (i.e., as those excited by 
illuminating with light incident perpendicularly to the 
graphene). But now, we show in Fig. |3^c),(d) the full 
plasmon dispersion relation as a function of frequency uj 
and wave vector fey parallel to the ribbon, both for inho- 
mogeneous doping produced by an external uniform field 
[Fig. 3'c)] and for a ribbon with uniform doping density 
[Fig. 3'd)]. The dispersion relations are rather different 
in both situations, with the inhomogeneous ribbon show- 
ing a denser set of modes, as well as more localization in 
the lowest-energy plasmons for large fey, as we show in 
Appendix [C] by means of near- field plots for the lowest- 
energy modes of both types of ribbons. 

Finally, let us mention that the inelastic plasmon de- 
cay rate is given by 7 within the Drude model in uni- 
formly doped structures However, 7 depends on po- 
sition for inhomogeneous doping. Using the DC mo- 
bility /i, one can estimate 7 = ev F /fi \Ef(x)\ ^ 2 x 
10 12 s _1 for Ep = 0.5 e V an d a typical measured mobil- 
ity fi = 10, 000 cm 2 /VsP^EI] Noticing that the local con- 
tribution to inelastic losses is proportional to Re{a} « 
(e 2 / , 7rti 2 )\EF\'y I 'uj 2 (i.e., independent of x), we conclude 
that the inhomogeneity of 7 is however translated into a 
uniform spatial distribution of losses. 

Conclusions.- We have shown that the plasmons of 
doped graphene ribbons are highly sensitive to the in- 
homogeneities of the doping charge density produced by 
realistic electrostatic landscapes. The doping profile can 
be engineered by adjusting the configuration of the gates 
relative to the graphene. We find an interesting scenario 
when a uniform external electric field is used to dope the 
graphene, leading to plasmons with very different char- 
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acteristics (e.g., induced charges piling up near the cen- 
ter of the ribbon) compared to those of uniformly doped 
graphene (in which plasmons pile up at the edges). This 
configuration can be used to avert losses associated with 
nonlocal effects at the edges, which are expected to be 
significant P=l The present study can be straightforwardly 
extended to other configurations, such as finite graphene 
nanoislands exposed to either backgates or side gates. 
Electrostatic charge accumulation at sharp edges can of- 
fer an additional handle to manipulate plasmon modes. 
In addition to the possibilities explored in this paper, one 
can use biased tips to produce localized disk-like doping 
areas at designated positions targeted by simply moving 
the tips above a graphene flake. In conclusion, the de- 
sign of electrostatic landscapes becomes a useful tool to 
engineer graphene plasmons. 
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ish MICINN (MAT2010-14885 and Consolider Nano- 
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4-248909-LIMA and FP7-ICT-2009-4-248855-N4E). 



Appendix A: Electrostatic calculation of doping 
surface charge distributions 



We consider three different geometries for electrically 
doping graphene, as illustrated in Fig. [4j Geometries 
in Figs. [4ja),(c) refer to doping of biased graphene rib- 
bons, which accumulate a net electric charge. In con- 
trast, the ribbon of Fig. [4^b), exposed to an external 
uniform field, has zero total charge. The charge distribu- 
tion in a backgated ribbon [Fig. Qa)] has been reported 
in the past p™^ particularly in the context of microstrip 
technology] 14 1 15 1 but we provide here a general procedure 
to calculate it based upon boundary elements, which we 
extend to other geometries under consideration. 
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1. Backgated graphene ribbon 

A ribbon placed at a potential V relative to a backgate [Fig. Qa)] displays a charge density —en(x) that can be 
calculated using the method of images (notice that n is the doping electron density). The charge depends on the 
coordinate across the ribbon x, which varies in the range < x < W where W is the width. Using the method of 
images, this problem is equivalent to two parallel ribbons vertically separated by a distance 2d and placed at potentials 
V (upper ribbon) and —V (lower ribbon), so that the backgate plane (z = 0) is at zero potential. The lower ribbon 
is thus represented by a charge density en(x). The potential at x in the upper ribbon is then given by 



V 



w 



dx I dy [— en{x')\ 



yj(x - x') 2 + y 2 y/(x - x') 2 + y 2 + 4d 2 



(Al) 



Analytically performing the integral along the y coordi- 
nate (perpendicular to the plane of Fig. [4| and using the 
notation r] = W/d, = x/d, and 

u = —V/ed, 



Eq. ( Al ) reduces to 

u= P d6 , n(0 , d)F(6,0'), 
Jo 



(A2) 



where 



F(0,0') =ln 



(0 - 0') 2 



We solve this integral equation by discretizing through 
a set of N equally spaced points 6j = (j + 1/2)7/ /iV, with 
j = 0, . . . , N — 1. Equation ( A2) is then approximated 
as 



U ~Y1 n (°j /d ) M n'i 



(A3) 



where 



r 



+r)/2N 



rj/2N 



d0' F(0j,0') 



(A4) 



is an integral over the interval surrounding point . We 
have assumed that n(x) is a smooth function, and al- 
though we find later that it diverges as ex 1/ yjx near the 
edges (see main paper), this divergence is integrable and 
contributes negligibly to the total integral for A 7 " ^> 1. 
From here, the charge distribution is found by inverting 
the matrix M , so that 



n(0jd) =u^2[M~ 



*33 



(A5) 



In practice, this method converges for N ~ 100. Similar 
conv ergen ce is obt ained for the geometries considered in 
Sees. A 2 and A3 Each curve in Fig. 1(a) of the main 
paper is actually consisting of two curves obtained with 
N = 100 and N = 500, and one cannot tell the difference 
between them on the scale of the plot. Finally, notice 
that the uniform electron density in the W/d ^> 1 limit 
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FIG. 4: Geometries under consideration for electrostatic doping of graphene: (a) single ribbon placed at a potential V relative 
to a backgate; (b) single ribbon exposed to a uniform external electric field; (c) two co-planar parallel ribbons set at opposite 
potentials, —V and V. 




x/W 

FIG. 5: Electrostatic doping electron density in backgated 
ribbons under the configuration sketched in Fig. [4|a) for dif- 
ferent values of the ratio of the ribbon width W to the gap 
distance d. The electron density n is normalized to the value 
n°° = —V/ Aired corresponding to the W/d ^> 1 limit. 
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FIG. 6: Electrostatic doping electron density across pairs of 
ribbons for different ratios of the ribbon width W to the gap 
distance d. The ribbons are placed at potentials —V and V, 
respectively (see inset), and the electron density n is normal- 
ized to the value n°° = —V/Aned. 



is given by n°° = u/Air. Also, we obtain the average 
of the Fermi energy (\Ep\) normalized to the W/d ^> 1 
limit Ef hv F s/\u\/4 as 



1\Ef\) 



3 



4n\n(8jd)/u\ 



where vp ~ 10 6 m/s is the Fermi velocity. (Inciden- 
tally, we consider the absolute value of Ep because the 
graphene response is nearly insensitive to the sign of Ep.) 

Figure [5] shows examples of doping charge densities 
from which we have extracted the Fermi energies shown 
in the inset of Fig. 1(a) of the main paper. 



Two co-planar parallel ribbons set at opposite 
potentials 



We can repeat the same analysis as in Sec. A 1 for 
two ribbons arranged as shown in Fig. [4jc) and set at 
potentials —V (left ribbon) and V (right ribbon). The 
separation between ribbons is d, and x = is chos e n at 
the left edge of the right ribbon. Equations (A2)-(A5) 



remain valid when we consider the charge d istri bution in 
the right ribbon, but now the kernel of Eq. ( A2) becomes 



F(0,0') = 21n 
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FIG. 7: Average Fermi energy (\Ef\) as a function of d/W 
for pairs of ribbons under the configuration of Fig. [4^ c). The 
Fermi energy is normalized to Ep° — Hvf\/\V\/ '4ecL 




FIG. 8: Doping charge density in a graphene ribbon subject 
to an external uniform electric field S. 



The charge density in the left ribbon is found from the 
symmetry n(x) = —n(—x — d). 

Figure [6] shows examples of doping charge densities 
from which we have extracted the Fermi energies shown 
in Fig. 2(a) of the main paper. The average Fermi energy 
(\Ef\) is shown in Fig. [7] as a function of d/W. 



3. Graphene ribbon in a uniform electric field 

For a ribbon subject to an external uniform electric 
field, only the component parallel to the ribbon £ can 
produce charge redistributions and local doping [see Fig. 
|4^b)]. Following the same procedure as in Sec. [XT] and 
taking the ribbon to be placed at zero potential, we can 



/W/2 noo 
dx' / dy[-en(x')] 
-W/2 J-oo 



y/(x- x') 2 + y 2 ' 



where the first term is the scalar potential produced by 
the external field. Now, using the normalization = 
x/W, the above equation reduces to 

,1/2 

£0= / d9'[-en(6'W)]F(6,6'), 

J -1/2 



where 



F(0,0') = ln 



(0 - Q') 2 



We use the discretization Oj = -1/2 + (j + 1/2) /N to 



write an expressio n si milar to Eq. ([A3]), where Mjj> is 
still given by Eq. (A4) with 77 = 1. Finally, the electron 
density is obtained from 

n(e j d) = (-£/e)J2[M- 1 ] jj ,e r . 

The doping charge density obtained from this equation 
is represented in Fig. [8] From here, we find a local Fermi 
energy Ep(x) = frvp \/7i\n{x)\ as shown in Fig. 3(a) of 
the main paper. The average Fermi energy is found to 
be (\E F \) = Q.6hv F y/\£\/e. 



Appendix B: Scaling of plasmon frequencies in the 
Drude model 



The electric scalar potential <p associated with the plas- 
mon modes in a doped graphene planar nanostructure 
satisfies the self-consistent electrostatic equation 



(R) 



d 2 R' 



IR-R 



- V • [<t(R» V0(R')] , (Bl) 



where an exp(— iuot) time dependence is undertood and 
all quantities are evaluated at the graphene plane (i.e., 
coordinates R and R/ are in that plane). This equa- 
tion represents the field induced by the self-consistent 
surface-charge density, which is calculated from the con- 
tinuity equation as (— i/uo)V • j in terms of the current j, 
and this is in turn proportional to the in-plane electric 
field via — <rV0, where a is the conductivity. Incidentally, 
a = outside the graphene, so that a presents a jump 
at the graphene edge, the grad ient of which generates an 
edge charge density. Equation ([Bl]) is valid for small pat- 
terning lengths compared to the light wavelength. The 
constant x takes the values x = 1 for suspended graphene 
in vacuum, x = 1/e for graphene embedded in a uniform 
dielectric of permittivity e, or x = 2/(e + 1) for graphene 
supported on a substrate.^ Now, we consider the Drude 



7 



model for the conductivity, 



<j(R, uj) 



e 2 i\E F (R)\ 
irh 2 (uj + ij) 



(B2) 



where < 7 <C uj and the R dependence comes from the 
inhomogeneous electron density n(R) (see Sec.[A]), which 
locally situates the Dirac points at an energy Ep(H) = 
Hvp \/7r|n(R)| relative to the Fermi levelP^ 

From the analysis of Sec. [Aj we can write the Fermi 
energy distribution as Ep(H) = (\Ep\) /(R), where / is a 
dimensionless envelope function. Likewise, distances can 
be scaled with a characteristic length W as U = R/W. 
Then, Eq. (Bl) becomes 



4>{\J) = c 



d 2 U' 
IU-U' 



■Vu-[/(U')Vu/^(U')]. 



where 



c 



-e 2 X 



7rh 2 W uj(uj + 27) 



is a dimensionless eigenvalue. The modes of the system 
satisfy this equation for specific choices of £, and there- 
fore, we conclude that the plasmon energies co> can be 
naturally normalized to the frequency 



oj = (e/hWx(\E F \)/W. 



(B3) 



Obviously, uj/ujq does not depend on the specific choice 
of (\Ep\) (the doping level), W (the size of the system), 
or x (the dielectric environment), and consequently, we 
present results normalized in this way in the main pa- 
per, which are universal for the kind of geometries under 
consideration, provided we stay within the limits of va- 
lidity of the Drude model (i.e., uj < Ep, and Ep smaller 
than the optical phonon energy ~ 0.2 eV^J). Addition- 
ally, because the electrostatic eigensystem is Hermitianp^ 



the eigenvalues ( are real, and therefore, the plasmon fre- 
quencies have imaginary part Im{co>} = —7/2, so that the 
plasmon lifetime is I/7, independent of geometrical and 
physical parameters within the Drude approximation. 



Appendix C: Computation of plasmon frequencies 
and near fields 



We solve the electrostatic problem of Eq. (Bl) by 
describing the graphene as a periodic array of surface 
dipoles with a small period compared to the characteris- 
tic lengths of the structure. This is the discrete surface- 
dipole approximation^ (DSDA), in which the polariz- 
ability of each element is taken such that a layer formed 
by a uniform lattice of dipoles has the same conductiv- 
ity as a uniform layer of graphene. The sum over dipole 
elements along y is performed before a self-consistent so- 
lution is sought, and therefore, the numerical problem 
reduces to solving a set of 2N linear equations with 2N 
variables (the dipole components along both x and y di- 
rections), where N is the number of dipoles across x. 
In practice, convergence is achieved with a few hundred 
dipoles for the dimensions considered in this work. Here, 
we have modified this method by allowing each element 
to depend through Ep on the spatial position along x. 

As an example of calculation, we show in Fig. |9ja),(b) 
plasmon dispersion diagrams for both an inhomoge- 
neously doped ribbon under the same conditions as in 
Fig. [8] and a homogeneously doped ribbon. These plots 
are the same as Fig. 3(c),(d) of the main paper. We 
also show the near fields calculated for the three lowest- 
energy modes in each case, with k\\W = 10. These are 
plasmons of monopole (c,f), dipole (d,g), and quadrupole 
(e,h) character. In the inhomogeneous ribbon, the modes 
are associated with large field enhancement near the cen- 
ter. In contrast, the uniform ribbon hosts modes with 
large field enhancement near the edges. 
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